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Abstract 

Heavy fourth generation Majorana neutrino can be stable and contribute to a 
small fraction of the relic density of dark matter (DM) in the Universe. Due to its 
strong coupling to the standard model particles, it can be probed by the current 
direct and indirect DM detection experiments even it is a subdominant component 
of the whole halo DM. Assuming that it contributes to the same fraction of the 
local halo DM density as that of the DM relic density in the Universe, we show that 
the current XenonlOO data constrain the mass of the stable Majorana neutrino to 
be greater than the mass of the top quark. In the mass range from 200 GeV to 
a few hundred GeV, the effective spin-independent cross section for the neutrino 
elastic scattering off nucleon is insensitive to the neutrino mass and mixing, and 
is predicted to be ~ 1.5 x 10~^^ cm^, which can be reached by the direct DM 
detection experiments soon. In the same mass region the predicted effective spin- 
dependent cross section for the heavy neutrino scattering off proton is in the range 
of 2 X 10-"^° cm^ ~ 2 X 10"^^ cm , which is within the reach of the ongoing DM 
indirect search experiments such as IceCube. We demonstrate such properties of 
the heavy neutrino DM in a fourth generation model with the stability of the fourth 
Majorana neutrino protected by an additional generation-dependent U{1) gauge 
symmetry. 
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1 Introduction 



Models with chiral fourth generation fermions are among the simplest and well-motivated 
extensions of the standard model (SM) and have been extensively studied [1]. The 
condition for CP symmetry violation in the SM requires at least three generations of 
fermions [2]. However, there is no upper limit on the number of generations from the 
first principle. In the SM the amount of CP violation is not large enough to explain 
the baryon-antibaryon asymmetry in the Universe. The inclusion of fourth generation 
quarks leads to two extra CP phases in quark sector and possible larger CP violation [3] , 
which is helpful for electroweak baryogenesis. With very massive quarks in the fourth 
generation, it has been proposed that the electroweak symmetry breaking may become 
a dynamical feature of the SM [1H7]. 

The recent LHC and Tevatron experiments have imposed constraints on the mass 
of the fourth generation quarks from possible pair production processes. For instance, 
the lower limit on the mass of the fourth generation up-type quark is found to be 
> 450 GeV from the search for the process ^4^4 — )■ WbWb — j- h?>jl^$j,, and the 
limit on the mass of the fourth generation down-type quark is m^^ > 490 GeV from 
the process ^4^4 — )■ WtWt — )■ £^£^h?>i$'j' [8]. The direct searches for extra quarks at 
Tevatron have set lower limits on the masses of U4 to be m„4 > 335 GeV, with the 
assumption that the mass splitting between M4 and the fourth generation down quark ^4 
is smaller than the W mass and the branching ratio of M4 — ?■ Wq is 100% [9]. The lower 
limit for the mass of is 771^^ > 385 GeV from the search for d^ — )■ Wt and the search 
for AW final states from the pair production of pp d^d^ [10], with the assumption 
that m„4 > rridA. Note that these limits are obtained with the assumptions of maximum 
mixing between the fourth and the third generation quarks and 100% branching ratio of 
the decay processes. The limits can be significantly weaker in models with suppressed 
mixing and decay branching ratio. 

In fourth generation models, the production of the Higgs boson can be enhanced by 
a factor of 5-9 due to the presence of two additional fourth generation quarks in the 
one- loop gluon-gluon fusion process. In a combined analysis of ATLAS and CMS, the 
mass of the Higgs boson has been excluded in the range 120-600 GeV which is based 
on the searches for — )■ 77, and Z^Z^{Z*^) in fourth generation models 

[TT] . In obtaining such a bound, it is assumed that the fourth generation fermions are 
heavy and do not contribute to the total width of the Higgs boson. Note that for a 
Higgs boson heavier than 600 GeV, the self-interaction of the Higgs boson may be very 
strong and even nonperturbative, while a Higgs boson lighter than 120 GeV may cause 
the problem of vacuum instability if the SM with fourth generation is valid up to the 
Planck scale. 
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The constraints on the masses of the fourth generation leptons are much weaker. The 
current lower bound on the mass of the unstable fourth generation charged lepton 64 is 

> 100.8 GeV from the search for 64 decaying into the fourth generation neutrino 1^4 
and boson. From the invisible width of boson, the lower bounds for the mass of an 
unstable is set to be 171^4^ > (101.3, 101.5, 90.3) GeV from the decay 1^4 — )■ (e, /i, r)iy^ 
in the case that is of Dirac type. If z/4 is of Majorana type the corresponding bounds 
are modified to be > (89.5,90.7,80.5) GeV. The constraints are much weaker if 
1^4 is a long-lived or stable particle. In this case the lower bound is around ~ mz/l: 

> 45.0(39.5) GeV for Dirac (Majorana) neutrino [I2]. A long-lived fourth neutrino 
can also relax the constraints from the precision electroweak data [13]. If the fourth 
generation neutrino is stable, it can be a potential candidate for the dark matter (DM) 
in the Universe. 

Heavy stable neutrinos with mass greater than ~ 1 GeV are possible candidates for 
the cold DM p31[T5] . However, if the neutrino is the dominant component of the halo DM, 
the current DM direct search experiments have imposed strong constraints on its mass. 
For instance, it has been shown that the spin-independent (SI) cross sections for the 
stable Dirac neutrino elastic scattering off nucleus can be a few order of magnitudes larger 
than the current DM direct search upper bounds due to the and exchanges [T6] - 
[T9] . In the case of Majorana neutrino, the spin-dependent (SD) cross sections for the 
elastic scattering can be very large. Previous analysis based on the assumption that the 
Majorana neutrino dominates the local halo DM density have ruled out the mass range 
10 GeV — 2 TeV from the upper bounds on the SD cross section [20]. The mass of the 
Majorana neutrino is also constrained from SI scattering through /i°-exchange if it has 
both Dirac and Majorana mass terms in the flavor basis [21]. On the other hand, it is 
well-known that for a neutrino heavier than ~ m^/2, the cross section for its annihilation 
is in general too large to reproduce the observed DM relic density. If the heavy neutrino 
mass is in the range m^/2 < 771^4^ < myi/, the annihilation into light fermion pairs // 
through s-channel Z° exchange contributes to a very large cross section. For the neutrino 
heavier than mw, the contribution from // channels decrease rapidly. However, other 
channels such as VF^IV^, Z^h^ etc. are opened. For these processes the corresponding 
cross section does not decrease with the increasing of the neutrino mass, resulting in 
a relic density always decreases with the growing of the neutrino mass, and a thermal 
relic density far below the observed total DM relic density [22]. Since the neutrino DM 
can only contribute to a small fraction of the relic density of DM, it is expected that it 
contributes to a small fraction of the halo DM density as well, and the two fractions are of 
the same order of magnitude. From model building point of view, it is easy to construct 
multi-component DM models with the heavy neutrino being a subdominant component 
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in the halo. Despite its very low number density in the halo, it can still be probed by the 
underground DM direct detection experiments due to its strong coupling to the target 
nuclei, which provides a way to search for new physics beyond the SM complementary 
to the LHC. In this case, the event rate of the DM-nucleus elastic scattering will depend 
on both the relic density and the cross section for the elastic scattering process. Since 
both of them have nontrivial dependence on the neutrino mass, the above mentioned 
constraints on the neutrino DM could be modified significantly. 

In this work, we explore the consequence of this possibility in a model with a fourth 
generation Majorana neutrino DM. The stability of the fourth Majorana neutrino pro- 
tected by an additional generation-dependent U{1) gauge symmetry. In the model the 
gauge-anomalies generated by the first three generation fermions are canceled by the ones 
from the fourth generation. We perform an updated analysis of the Majorana neutrino 
DM in light of the recent Higgs search results at the LHC and the recent DM direct de- 
tection results such as XenonlOO and SIMPLE etc.. The relic density of the neutrino DM 
is obtained by calculating the annihilation cross sections for all the possible final states 
such as //, W^W"^, Z^Z'^, Z^h^ and h^h'^ etc.. Under the assumption that the neutrino 
contributes to the same fraction of the local halo DM density as that of DM relic density 
in the Universe, we calculate the effective DM-nucleus elastic scattering cross sections 
which are the cross sections rescaled by the fraction of the halo DM density contributed 
by the heavy Majorana neutrino. The results show that the current XenonlOO data 
constrain the mass of the Majorana neutrino to be greater than ~ 175 GeV. In the 
mass range from 200 GeV to a few hundred GeV, the spin-independent cross section for 
the neutrino elastic scattering off nucleon is predicted to be ~ 1.5 x 10~^^ cm^, which is 
insensitive to the neutrino mass and can be reached by the direct DM search experiments 
in the near future. In the same mass range the predicted spin-dependent cross section 
for neutrino proton scattering is in the range 2 x 10"^'' cm^ ~ 2 x 10~^^ cm^, which is 
within the reach of the ongoing IceCube experiment. Although the analysis is performed 
in a specific model, the results are valid for a wide range of models with heavy Majorana 
neutrino DM. 

This paper is organized as follows, in section [21 we present the details of the model 
with stable fourth generation Majorana neutrino and the interactions relevant to the 
calculation of the relic density and DM-nucleus elastic scattering cross sections. In 
section [3|, we calculate the relic density of Majorana neutrino DM. In section |U we 
give the prediction for spin-independent and spin-dependent cross sections for the heavy 
neutrino elastic scattering off nucleon and compared them with the latest DM search 
experiments. Finally, we give the conclusions in section |5l 
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2 A model with stable fourth generation neutrino 



In this work we consider a simple extension of the SM with a sequential fourth generation 
and an additional U{1)f gauge symmetry. The U{1) extensions to the SM are well 
motivated from the point view of grand unification such as the 5*0(10) and Eq and have 
rich phenomenology [23] which can be reached by the on going LHC experiments. The 
flavor contents in the model are given by 

QiL = {^'^ ] ^ ^iL = { I , UiR, diR, UiR, em (z = 1, . . . ,4). (1) 
\diL J \eiL J 

All the fermions in the model are vector-like under the extra gauge interactions associated 
with U{1)f- The U{1)f charges of the fermions could be generation- dependent. In order 
to evade the stringent constraints from the tree-level flavor changing neutral currents 
(FCNCs), the U{1)f charges Qgi for the first three generation quarks are set to be the 
same, i.e. Qgi = Qg, {i = 1, 2, 3) while Qg^ = —3Qg for the fourth generation quarks. 
Similarly, the U{1)f charges for the first three generation and the fourth generation 
leptons are Ql and —3Ql, respectively. In general, Qg and Ql can be different. For 
simplicity, in this work we take Qg = Ql = 1. With this set of flavor contents and U{1)f 
charge assignments, it is straight forward to see that the new gauge interactions are 
anomaly-free. Since the gauge interaction of U{1)f is vector-like, the triangle anomalies 
of [U{l)Ff, [SU{3)c]'^U{1)f and [gravity] ^[/(l)^ are all vanishing. The anomaly of 
U{1)y[U{1)f]'^ is zero because the f/(l)y hypercharges cancel for quarks and leptons 
separately in each generation, namely J2i~'^qL + ^qr) = and X^l + X^i?) = 0. The 
anomaly of [SU{2)l]'^U{1)f is also zero due to the relation 

4 4 

Y^Qgi = and Y,QL^ = 0■ (2) 

1=1 1=1 

Thus in this model, the gauge anomalies generated by the first three generation fermions 
are canceled by that of the fourth generation one, which also gives a motivation for the 
inclusion of the fourth generation. 

The gauge symmetry U{1)f is to be spontaneously broken by the Higgs mechanism. 
For this purpose we introduce two SM singlet scalar fields (j)a,b which carry the U{1)f 
charges Qa = —^Ql and Qb = 6Ql respectively. The U{1)f charges of (f)a,b are arranged 
such that (pa can have Majorana type of Yukawa couplings to the right-handed neutrinos 
of the first three generations (z = 1, 2, 3) while (pb only couples to the fourth generation 
neutrino z/4/j. After the spontaneous symmetry breaking, the two scalar fields obtain 
vacuum expectation values (VEVs) {4>a,b) = Va,b/V^- 
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The relevant interactions in the model are given by 



1 

' l_K7- J \-TJ -7-7-/ 2 



-Yl^vtR<t^am = 1, 2, 3) - -Yriy^R'Pb'^m - Vi<Pa, 06, H) + H.c. (3) 



where fi stand for left- and right-handed fermions, and H is the SM Higgs doublet. 
Dfj,fi = {d^ — igiT°'WI^ — iYg2Bn — iQ jgpZ'^fi is the covariant derivative with Z'^ the 
extra gauge boson associated with the U{1)f gauge symmetry, and gp the corresponding 
gauge coupling constant. Since (f)a,b are SM singlets, they do not play any role in the 
electroweak symmetry breaking. Thus Z' obtains mass only from the VEVs of the scalars 

ml = 9l{Qlvl + Qlvll (4) 

and it does not mix with Z^ boson in mass term at tree level. It can only mix with Z^ 
through kinetic mixing. In this work we assume that the effect of the kinetic mixing is 
small and negligible. 

From the U{1)f charge assignments in the model, the four by four Yukawa coupling 
matrix is constrained to be of the block diagonal form 3® 1 in the generation space. Since 
the U{1)f charges are the same for the fermions in the first three generation, there is no 
tree level FCNC induced by the Z'-exchange in the physical basis after diagonalization. 
Thus a number of constraints from the low energy flavor physics such as the neutral 
meson mixings and the 6 — )■ 57 can be avoided. 

The contribution to the muon g — 2 from the Z' boson at one-loop can be estimated 
as pi] 

The current experimental data requires that Aa^ < 3.9 x 10"^ [2H1E6], which can be 
translated into a lower bound on the VEVs of 0a,fe: V Qlvl + Ql'^l — 1-45 x 10^ GeV. 
Note that the dependence on the coupling constant (^f is canceled by the one in the mass 
of Z' in the expression of Aa^. The direct search for the process e~^e~ — — t- at 
the LEP-II leads to a lower bound on the ratio of the mass to the coupling to leptons: 
Mz'/gp > 6 TeV [27] for vector-like interactions, which corresponds to a more stringent 
lower bound: ^/QlvlTQM > 6 TeV. 

The current searches for narrow resonances in the Drell-Yan process pp ^ Z' ^ 
at the LHC impose an alternative bound on the mass and the couplings of the Z' boson. 
In the narrow width approximation, the cross section for the Drell-Yan process can be 
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parametrized as a = {n/ {4:8s))[cuWu{s, M"^,) + CrfWd(s, M|,)], where s is the squared 
center of mass energy and Wu,d{s, M^,) are model-independent functions depending on 
Mz' [2Z]- The coefficients Cu,d are related to the Z' couplings to the quarks as Cq = 
+ Q^FdBri^Z' — )■ (g = u,d) where QqL,qR are the charges of the left- and 

right-handed quarks of a generic U{1) gauge symmetry with coupling strength g. The 
limit on the cross section reported by the experiments can be recast model-independently 
as limit contours in the Cu — Cd plan with unique contour for each value of Mz' |27]. Using 
the Cu — Cd contours, the limits on one type of Z' model can be translated into that of 
other models. For a model with sequential neutral gauge boson Zg^^j which by definition 
has the same couplings as that for the SM Z^ boson [23], the latest lower bounds on its 
mass Mz'^^^^ is 1.94 TeV from CMS [28] and 1.83 TeV from ATLAS [29] respectively. In 



the Z'ss;^^ model cf^^^ ^ 2.36 x 10"^ and cf^*^ ^ 3.66 x 10"^ The bound on M^^^^^ can 
be translated into the bound on the mass and couplings of the Z' in this model in which 
c„ = Q = lSglBr{Z' /i+yU"). By requiring that c„ < c^^^' at Mz' = 1.94 TeV, we 
obtain gp < 0.051. Numerical analyses have shown that for a given limit on the cross 
section, as the combination CuWu + CdWd increases or decreases by an order of magnitude, 
the mass limit changes roughly by 500 GeV [30l[3T]. Thus as a rough estimation one can 
obtain limits on gp for other values of Mz', for instance gp < 0.029 for Mz' = 1.44 TeV 
and gp < 0.0051 for Mz' = 0.94 TeV, respectively. 

In this model, there is no mixing between the fourth and the first three generation 
quarks, thus the null results in searching for U4 — > bW^ and — )■ tW^ at the LHC do 
not impose constraints on the masses of M4 and d^. The model is not constrained by the 
FCNC processes U4 — j- tX or d^ — )■ bX either [32] . 

The fourth generation neutrinos obtain both Dirac and Majorana mass terms through 
the vacuum expectation values (VEVs) of H and (pb- In the basis of (i^l, i^'f>)'^ the mass 
matrix for the fourth neutrino is given by 



{ , (6) 




where rriD = YI'vh/V^ with vh = 246 GeV and niM = yr"^<i>bl The left-handed 
components {vyp v^^l) '^^^ ^"^0 mass eigenstates are related to the ones in the flavor 
eigenstates by a rotation angle B 



IL 



-i[cevp - sgu"^), v^^' = seVL + cgz/^, (7) 



where sg = sin 6 and cg = cos 6'. The value of 6 is defined in the range (0,7r/4) and is 
determined by 

tan2e = ^, (8) 
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with ^ = (vr/4) corresponding to the hmit of minimal (maximal) mixing. The phase i 
is introduced to render the two mass eigenvalues real and positive. The two Majorana 
mass eigenstates are xi = i^il^ + '^il^^ and X2 = ^^l + ^^l" ^ respectively. The masses of 
the two neutrinos are given by m\^2 = ( + 4m|, =]= mA/)/2. In terms of the mixing 
angle 6 they can be rewritten as 

with mi < 7712. Note that for all the possible values of 6 the lighter neutrino mass 
eigenstate xi consists of more left-handed neutrino than the right-handed one, which 
means that xi always has sizable coupling to the SM boson. Therefore the LEP-II 
bound on the mass of stable neutrino is always valid for xi, which is insensitive to the 
mixing angle. 

As the fermions in the first three generations and the fourth generation have different 
U{1)f charges, the fourth generation fermions cannot mix with the ones in the first 
three generations through Yukawa interactions. After the spontaneous breaking down 
of U{1)f, there exists a residual Z2 symmetry for the fourth generation fermions which 
protect the fourth neutrino Xi to be a stable particle if it is lighter than the fourth 
generation charged lepton 64, which makes it a possible dark matter candidate. 

Although the stability of xi can be achieved easily by imposing an ad hoc Z2 sym- 
metry only on the fourth generation, it is theoretically more attractive to consider a Z2 
which originates from a U{1) gauge symmetry in which all the fermions are charged. 
First, models with extra U{1) gauge symmetries are well motivated from various GUT 
models such as S'O(IO) and Eq etc., and it is common to link a discrete Zj^ symmetry 
to a broken continuous U{1) symmetry in model building. Second, in this model all the 
four generations are treated in parallel. The only difference is the U{1) charge of the 
fourth generation fermions. Due to this U{1) charge difference, after the spontaneous 
symmetry breaking by the VEVs of (pa and (pb , the U{1) is broken into two separate 
Z2 symmetries, with one on the first three generation and the other one on the fourth 
generation which stabilizes the fourth generation neutrino. Third, in this model setup, 
all the fermions are vector-like under the extra U (1) gauge interactions. The cancellation 
of gauge anomalies automatically requires the existence of the right-handed neutrinos in 
each generation. Thus in this model it is natural to have extra stable neutrinos. 

In the mass basis the interaction between the massive neutrinos and the SM Z^ boson 
is given by 

Cnc = , ^\ [-clxil^l'xi - SeX27'^7'x2 + 1iceSeXil''X2] Z,, (10) 

4 cos (7w 
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where gi is the weak gauge couphng and 6w is the Weinberg angle. The Yukawa inter- 
action between Xi.2 and the SM Higgs boson is given by 

£y = -H^ I ^ I [cgSeXiXi + ceSeX2X2 - i^c] - sl)xil^X2\ ■ (H) 
vh \se J 

Since all the fermions are vector-like under the U{1)f gauge interaction and Xi,2 are 
neutral particles, only the off-diagonal interaction X1X2Z' is allowed. 

3 Annihilation cross sections and relic density of the 
fourth generation neutrino dark matter 

The thermal relic density of Xi is related to its annihilation cross section at freeze out. 
The Feynman diagrams for all possible annihilation channels are shown in Fig. [TJ When 
the mass of xi is smaller than that of the boson, XiXi can only annihilate into 
light SM fermion pairs through s-channel /h^ exchange. For Majorana neutrino the 
annihilation cross sections are suppressed by the small masses of the final state fermions. 
However, large enhancement of the annihilation cross section occurs if the mass of xi 
is close to rtizjl such that the intermediate is nearly on shell. When the fourth 
generation neutrino is heavier than the boson, the ly^H^^ channel will open and 
become important. The process X\X\ ~^ W^W^ involves s-channel and /i^ boson 
exchange and t-channel 64 exchange. All these three intermediate states must be included 
in order to maintain the unitarity. For even heavier xi, the final states can be 
Z^hb^ and The Z^Z'^ channel contains the process of s-channel exchange and the 

t-channel Xi,2 exchange. Similarly, the Z'^hb^ channel contains s-channel Z° exchange and 
the t-channel xi,2 exchange processes. The channel is dominated by the t-channel 
processes due to the large Yukawa coupling. When x\ is heavier than the top quark, 
the tt final states must be included. The cross section can be large as the top quark has 
large Yukawa couplings to the Higgs boson. In this work we neglect the Z' Z' final states 
by assuming that Z' is always heavier than x\- 

We numerically calculate the cross sections for X\X\ annihilation into all the relevant 
final states using CalHEP 2.4 [33]. In order to determine the DM relic density, one needs 
to calculate the thermally averaged product of the DM annihilation cross section and 
the relative velocity 



where T is temperature and K\^2{p^^ are the modified Bessel function of the second kind. 




(12) 
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The relic abundance can be approximated by 



nn^ . '"J:"':' . (13) 



where x = rrii/T is the rescaled inverse temperature. Xi? ~ 25 corresponds to the 
decouphng temperature, g^, = 86.25 is the number of effective relativistic degree of 
freedom at the time of freeze out, and Mpi = 1.22 x 10^^ GeV is the Planck mass scale. 
In the left panel of Fig. [21 we show the value of (av) at the time of freeze out 
= 25 as function of mi. In the figure, the contributions from individual final states 
are also given. The peak at ~ rnz/2 corresponds to the case of resonant s-channel 
annihilation when the intermediate boson is nearly on shell. During the numerical 
calculations the mass difference between the charged fourth generation lepton 64 and xi 
is set to be rrie^ — m^-^ = 50 GeV. Note that a very large mass difference is subject to 
strong constraints from the electroweak precision data [M] . The results are found to 
be insensitive to the mass difference. In the left panel of Fig. [2] the mixing angle 6 is 
fixed at 30°. The recent LHC experiments have placed strong constraints on the mass of 
h^. For the fourth generation model the mass of is constrained to be below 120 GeV 
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or heavier than 600 GeV [TT]. In the numerical calculations the mass of the SM Higgs 
boson is fixed at 115 GeV. 

In the case that the Higgs boson is light enough to be among the final states of 
the XiXi annihilation, the contributions from the Z^h^ and h^h^ final states can be 
important. It has already been noticed that the Z^h^ channel can be as importance as 
the W^W^ channel in the case of heavy Dirac neutrino DM [19]. In the Majorana case, 
since the W^W^ cross section is strongly suppressed by the smallness of the relative 
velocity [35], the Z^h^ channel can give the dominant contribution to the cross section. 

In the right panel of Fig. [2], we show the quantity 

m - (14) 

which is the ratio of the relic density of xi to the observed total DM relic density 
^DAih"^ = 0.110 ± 0.006 [12] as function of the mass of Xi fo^^ different values of the 
mixing angle 6. The results show a significant dependence on the mixing angle 6. For 
smaller mixing angle 6 the couplings between Xi and gauge bosons W^, Z are stronger, 
resulting in a smaller relic density. The results also clearly show that due to the large 
annihilation cross section, xi cannot make up the whole DM in the Universe. Xi can 
contribute to ~ 20 — 40% of the total DM relic density when its mass is around 80 GeV. 
But for mi ^ nit, it can contribute only a few percent or less to the whole DM. 

However, since xi has strong couplings to and Z^, even in the case that the number 
density of xi is very low in the DM halo, it is still possible that it can be detected by 
its elastic scattering off nucleus in direct detection experiments. Given the difficulties 
in detecting such a neutral and stable particle at the LHC, there is a possibility that 
the stable fourth generation neutrino could be first seen at the DM direct detection 
experiments. 

In this model there exists three right-handed neutrinos Viji {i = 1,2,3) in the first 
three generations. In the physical basis there may exist three sterile neutrinos, provided 
that the mixings between the left- and right-handed neutrinos are tiny. If one of the 
sterile neutrinos has a mass around keV scale, it can be a good candidate for warm 
dark matter which could be the dominant component of the DM in the Universe. The 
warm dark matter may provide a solution to some of the known problems in the DM 
simulations based on cold DM, such as reducing the number of subhalos and smoothing 
the cusps in the DM halo center. In the SM with only right-handed neutrinos, the sterile 
neutrino may obtain the correct relic density through non-thermal production [36]. In 
this model, since the keV sterile neutrino can annihilate into the light active neutrinos 
through the extra U{1)f gauge interactions, it can also be a thermal relic. In this 
case, the correct relic density can be obtained by including the effect of entropy dilution 
through the decays of other heavier sterile neutrinos [37] . 
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Figure 2: Left) Thermally averaged product of the annihilation cross section and the 
relative velocity as function of the mass of the stable heavy neutrino Xi- The contribu- 
tions from individual final states are also shown; right) the rescaled Xi ^^^lic density tq 
as function of the mass of xi- The shaded region is excluded by the LEP-II experiments. 

4 Direct detections of the fourth generation neu- 
trino dark matter 

The generic formula for the differential event rate of DM-nucleus scattering per nucleus 
mass is given by 



where Eji is the recoil energy, is the scattering cross section corresponding to the zero 
momentum transfer, itidm is the mass of the DM particle, fiN = 'mDM'mN / (jnoM + "^Ar) 
is the DM-nucleus reduced mass, F{E^ is the form factor, and f{v) is the velocity 
distribution function of the halo DM. The local DM density puhi is often set to be equal 
to po — 0.3 GeV/cm"^ (for updated determinations of po see e.g. [38]) which is the local 
DM density inferred from astrophysics based on a smooth halo profile. Since the neutrino 
DM can only contribute to a small fraction of the relic density of DM, it is likely that 
it also contributes to a small fraction of the halo DM density, namely, its local density 
Pi is much smaller than po- If the DM particles are nearly collisionless and there is no 
long range interactions which are different for different DM components, the structure 
formation process should not change the relative abundances of the DM components. In 
this work, we assume that pi is proportional to the relic density of Xi iii the Universe, 




(15) 
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namely 



PO i 'DM 



or Tp ^ rn. Consequently, the expected event rates of the DM-nucleus elastic scattering 
will be scaled down by Vp. In order to directly compare the theoretical predictions with 
the reported experimental upper limits which are often obtained under the assumption 
that the local DM particle density is po? we shall calculate the rescaled elastic scattering 
cross section 

a = rpa^ r^a, (17) 

which corresponds to the event rate to be seen at the direct detection experiments. Note 
that a depends on the mass of Xi through the ratio even when a is mass-independent. 

The spin-independent DM-nucleon elastic scattering cross section in the limit of zero 
momentum transfer is given by [39] 

SI [ZU + {A-Z)Uf 

-j^ , (18) 

where Z and A — Z are the number of protons and neutrons within the target nucleus, 
respectively. /i„ = mim„/ (mi + m,„) is the DM-nucleon reduced mass. The coupling 
between DM and the proton (neutron) is given by 

hAn) - iTq — + T^JtG 2^ ^9^;^' (1^^ 

lltq ( IILq 

q=u,d,s ^ q=c,b,t ^ 

with /^^"'■' the DM coupling to light quarks and f^Q^ = 1 — ^q=u d s fj^q^ ■ ^^^^ 
that the elastic scattering is dominated by t-channel Higgs boson exchange, the relation 
/„ ~ fp holds and one has af^ ~ 4/^/i^/7r. In numerical calculations we take f^^ = 
0.020 ± 0.004, = 0.026 ± 0.005, = 0.118 ± 0.062, f^^ = 0.014 ± 0.003, /^^ = 
0.036 ± 0.008 and f^^ = 0.118 ± 0.062 [lO]. The coefficient aq in the model is given by 

S = (20) 

The value of aq is proportional to mi, thus larger elastic scattering cross section is 
expected for heavier xi- Note that in terms of mi the coefficient aq is proportional to Cg. 
Part of the mixing effects has been absorbed into the mass of xi- In the limit of 6' — 0, 
mi is approaching zero and the couping between xi and is vanishing as expected. The 
value of aq has a strong dependence on m/j. As the latest LHC data exclude the mass of 
in the range 120 GeV-600 GeV in the presence of fourth generation fermions [H], we 
fix rrih = 115 GeV in the numerical calculations. The quark mass rriq in the expression 
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of ttq cancels the one in the expression of fp{n)- Thus there is no quark mass dependence 
in the calculations. 

In Fig. [3] we give the predicted spin-independent effective cross sections a^' for 
the fourth generation neutrino elastic scattering off nucleon as function of its mass for 
different values of the mixing angle 6. One sees that even after the inclusion of the 
rescaling factor Vp, the current XenonlOO data can still rule out a stable fourth generation 
neutrino in the mass range 55 GeV < mi < 175 GeV which corresponds to < 1%. 
Thus the stable fourth generation neutrino must be heavier than the top quark, and can 
only contribute to a small fraction of the total DM relic density 

On the other hand, for m^^ > 200 GeV, the cross section does not decrease with 
increasing, and is nearly a constant a^^ ~ 1.5 x 10~'^^cm^ in the range 200 GeV < 
< 400 GeV. This is due to the enhanced Yukawa coupling between the fourth 
generation neutrino and the Higgs boson which is proportional to m^-^ , as it is shown in 
the expression of Og. One can see from the Fig. [3] that the result is not sensitive to the 
mixing angle 6 either, which is due to the compensation of the similar dependencies on 
6 in the relic density. For instance, the cross sections for the W^W^ and Z^Z^ channel 
of XX annihilation are proportional to Cg, which compensates the ^-dependence in the 
for the elastic scattering processes. 

The Majorana neutrino DM can contribute to spin-dependent elastic scattering cross 
section through axial-vector interaction induced by the exchange of the Z^ boson. At 
zero momentum transfer, the spin-dependent cross section has the following form [39] 

32 7 _|_ 1 

= -Glfil—— (a^iSp) + an{Sn)f , (21) 

TT J 

where J is the spin of the nucleus, ap(„) is the DM effective coupling to proton (neutron) 
and («S'p(„)) the expectation value of the spin content of the nucleon within the nucleus. 
Gp is the Fermi constant. The coupling ap(„) can be written as 

«.H = E7^^f^ (22) 

where dg is the DM coupling to quark and Ag*-"'' is the fraction of the proton (neutron) 
spin carried by a given quark q. The spin-dependent DM-nucleon elastic scattering cross 
section is given by 

- (tI;^""" ^ ^ tI;^- "') ^ ■ '''' 

In numerical calculations we take = 0.77, A^ = —0.40, A^ = —0.12 |13], and use the 
relations A" = A^, = A^, A" = A^. The coefficients dq in this model are given by 

Of 

du = -dd = -ds = (24) 
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Figure 3: Effective spin- independent cross section 5"f which is rescaled by Vp ^ 
for Xi elastically scattering off nucleon as function of the mass of Xi- Four curves 
correspond to the mixing angle 9 = 10°(sohd), 20° (dashed), 30° (dotted) and 40° (dot- 
dashed) respectively. The current upper limits from CDMS [H] and XenonlOO [42] 
experiments are also shown. 



For the axial-vector interactions, the coupling strengths do not depend on the electro- 
magnetic charges of the quarks. 

In Fig. IHwe show the predicted effective spin- dependent DM-neutron cross section 
^ as function of the neutrino mass for different mixing angles, together with various 
experimental upper limits. Since af^ is independent of m^-^, the dependency of a^^^-, 
on the neutrino mass comes from the dependency of on m^-^, which can be seen by 
comparing Fig. H] with Fig. |2l The XenonlO data is able to exclude the neutrino DM 
in the mass range 60 GeV < m^^ ^120 GeV, which is not as strong as that from the 
XenonlOO data on spin- independent elastic scattering cross section. For a heavy neutrino 
DM with mass in the range 200 GeV < m^^ < 400 GeV the predicted spin-dependent 
cross section is between 10~^° cm^ and 10~^^ cm^. 

In Fig. |5]we give the predicted spin-dependent DM-proton cross section dp^. The 
cross sections for Majorana neutrino DM scattering off proton and neutron are quite 
similar, which is due to the fact that the relative opposite signs in and are 
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compensated by the opposite signs in du and dn- So far the most stringent hmit on the 
DM-proton spin-dependent cross section is reported by the SIMPLE experiment [H]. 
The SIMPLE result is able to exclude the mass range 50 GeV < m^^ < 150 GeV, 
which is compatible with the constraints from XenonlOO. In Fig. |5l we also show 
the upper limits from indirect searches using up-going muons which are related to the 
annihilation of stable fourth generation neutrinos captured in the Sun. The limit from the 
Super-K experiment is obtained with the assumption that 80% of the DM annihilation 
products are from bb, 10% from cc and 10% from rr respectively |15]. In the range 
170 GeV < m^^ < 400 GeV, the limit from Super-K is ~ 5 x 10-^9 cm^. The IceCube 
sets a stronger limit a^^ < 2 x 10~^° cm^ for the DM mass at 250 GeV |46]. This limit 
is obtained with the assumption that the DM annihilation products are dominated by 
P^±iy=F jf f^i^Q annihilation products are dominated by bb, the limit is much weaker, for 
instance d-p^ < 5 x 10"^*^ cm^ for the DM mass at 500 GeV [16]. Note that in this model, 
the dominant final state is Z^h^. The expected limit should be somewhere in between. 
Nevertheless, the IceCube has the potential to test the predictions in this model. 

Different assumptions on the value of Vp and the nature of the heavy stable neutrino 
may result in different limits. For instance, in Ref. [20], an excluded mass range of 10 
GeV-2 TeV is obtained from the Xenon 10 data on the cross section of the spin-dependent 
DM-nucleus elastic scattering, which is based on the assumption that the local halo DM 
is entirely composed of stable Majorana neutrino, i.e. Vp = 1, and the neutrino has the 
same couplings to the boson as that of the SM active neutrinos. As in the present 
model we have Vp ^ tq <^ 1 and the coupling to the Z^ boson depends on the mixing 
angle, the resulting constraints are different significantly. 

5 Conclusions 

In conclusion, we have investigated the properties of stable fourth generation Majorana 
neutrinos as dark matter particles. Although they contribute to a small fraction of the 
whole DM in the Universe, they can still be easily probed by the current direct detection 
experiments due to their strong couplings to the SM particles. We have considered a 
fourth generation model with the stability of the fourth Majorana neutrino protected 
by an additional generation-dependent U{1) gauge symmetry. In the model the gauge- 
anomalies generated by the first three generation fermions are canceled by the ones from 
the fourth generation. We have shown that the current XenonlOO data constrain the 
mass of the stable Majorana neutrino to be greater than the mass of the top quark. For a 
stable Majorana neutrino heavier than the top quark, the effective spin-independent cross 
section for the elastic scattering off nucleon is found to be insensitive to the neutrino mass 
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Figure 4: Effective spin-dependent cross section 5"„ which is erf rescaled by ^ 
for Xi elastically scattering off neutron as function of the mass of Xi- Four curves corre- 
spond to the mixing angle 9 = 10°(sohd), 20° (dashed), 30° (dotted) and 40° (dot-dashed) 
respectively. The current upper limits from various experiments such as KIMS [47], 
CDMS [IE] and XenonlO [20] are also shown. 

and is predicted to be around 10"^^ cm^, which can be reached by the direct DM search 
experiments in the near future. The predicted effective spin-dependent cross section for 
the heavy neutrino scattering off proton is in the range 10~^° cm^ ~ 10^^^ cm^, which 
is within the reach of the ongoing DM indirect search experiments such as IceCube. 
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Figure 5: Effective spin-dependent cross section which is 0"^ rescaled by ^ tq 
for Xi elastically scattering off proton as function of the mass of Xi- Four curves corre- 
spond to the mixing angle 6 = 10°(sohd), 20°(dashed), 30°(dotted) and 40° (dot-dashed) 
respectively. The current upper limits from various experiments such as KIMS |47j . 
CDMS [IE], XenonlO [20], Coupp |19], Picasso [50], SIMPLE |5l], SuperK |15], and 
IceCube |16] are also shown. 
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